“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


2013-03 


WAVE AND CURRENT OBSERVATIONS INA 
TIDAL INLET USING GPS DRIFTER BUOYS 


Mcintyre, Stephen A. 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/32869 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


\§ D U DL EY research materials and institutional publications created by the NPS community. 
«iit Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed -- and published -- scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 





NAVAL 
POSTGRADUATE 
SCHOOL 


MONTEREY, CALIFORNIA 


THESIS 


WAVE AND CURRENT OBSERVATIONS IN A TIDAL 
INLET USING GPS DRIFTER BUOYS 


by 
Stephen A. McIntyre 


March 2013 


Thesis Advisor: Thomas H.C. Herbers 
Second Reader: Tim T. Janssen 





Approved for public release; distribution is unlimited 


THIS PAGE INTENTIONALLY LEFT BLANK 


Public reporting burden for this collection of information is estimated to average | hour per response, including the time for reviewing instruction, 
searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send 
comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to 
Washington headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 
22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
March 2013 Master’s Thesis 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 
WAVE AND CURRENT OBSERVATIONS IN A TIDAL INLET USING GPS 
DRIFTER BUOYS 


6. AUTHOR(S) Stephen A. McIntyre 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION 
Naval Postgraduate School REPORT NUMBER 
Monterey, CA 93943-5000 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
N/A AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official policy 
or position of the Department of Defense or the U.S. Government. IRB Protocol number N/A 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited 
13. ABSTRACT (maximum 200 words) 


The Wave Resolving Drifter (WRD) Buoy, developed at NPS, was refined by adding an accelerometer and utilizing 
horizontal Doppler velocity measurements to better resolve the wave surface motions. The WRD Buoy was validated 
against a Datawell Waverider-GPS buoy by comparing bulk wave statistics, wave energy spectra and wave directional 
spectra. The validation was performed in deep water, offshore of Monterey Bay, CA. Horizontal measurements from 
Doppler velocities showed improved capability, particularly in the wind-wave band. Vertical measurements were 
significantly improved through the addition of the accelerometer. A large array of WRD buoys was deployed in the 
Golden Gate channel at the entrance to San Francisco Bay and allowed to drift with the outgoing ebb current. The 
simultaneous deployment of many drifters provides a unique view of the details of the temporal and spatial evolution 
of the wave field as it propagated through regions of bathymetric and wave-current interactions. Wave heights 
increased as the wave field propagated over the ebb-tidal shoal and at the entrance to the channel, in the presence of a 
strong opposing ebb current, as expected. Inside the channel, strong dissipation of high frequency wind waves was 
observed in the opposing current. 


14, SUBJECT TERMS 15. NUMBER OF 
Ocean waves, Tidal currents, Wave buoys, GPS, Accelerometer, Ebb-tidal shoal, Wave-bathymetry PAGES 
interaction, Wave-current interaction 99 


16. PRICE CODE 


17. SECURITY 18. SECURITY 19. SECURITY 20. LIMITATION OF 
CLASSIFICATION OF CLASSIFICATION OF THIS CLASSIFICATION OF ABSTRACT 
REPORT PAGE ABSTRACT 
Unclassified Unclassified Unclassified UU 
NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 





THIS PAGE INTENTIONALLY LEFT BLANK 


i 


Approved for public release; distribution is unlimited 


WAVE AND CURRENT OBSERVATIONS IN A TIDAL INLET USING GPS 
DRIFTER BUOYS 


Stephen A. McIntyre 
Lieutenant Commander, United States Navy 
B.S., Auburn University, 2002 
M.S., University of Southern Mississippi, 2009 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN METEOROLOGY AND OCEANOGRAPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 


March 2013 
Author: Stephen A. McIntyre 
Approved by: Thomas H.C. Herbers 
Thesis Advisor 


Tim T. Janssen 
Second Reader 


Peter C. Chu 
Chair, Department of Oceanography 


ili 


THIS PAGE INTENTIONALLY LEFT BLANK 


iv 


ABSTRACT 


The Wave Resolving Drifter (WRD) Buoy, developed at NPS, was refined by adding an 
accelerometer and utilizing horizontal Doppler velocity measurements to better resolve 
the wave surface motions. The WRD Buoy was validated against a Datawell Waverider- 
GPS buoy by comparing bulk wave statistics, wave energy spectra and wave directional 
spectra. The validation was performed in deep water, offshore of Monterey Bay, CA. 
Horizontal measurements from Doppler velocities showed improved capability, 
particularly in the wind-wave band. Vertical measurements were significantly improved 
through the addition of the accelerometer. A large array of WRD buoys was deployed in 
the Golden Gate channel at the entrance to San Francisco Bay and allowed to drift with 
the outgoing ebb current. The simultaneous deployment of many drifters provides a 
unique view of the details of the temporal and spatial evolution of the wave field as it 
propagated through regions of bathymetric and wave-current interactions. Wave heights 
increased as the wave field propagated over the ebb-tidal shoal and at the entrance to the 
channel, in the presence of a strong opposing ebb current, as expected. Inside the channel, 


strong dissipation of high frequency wind waves was observed in the opposing current. 
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I. INTRODUCTION 


A. NEARSHORE WAVE PROCESSES 


Wind-generated surface waves are present on every body of water. In the open 
ocean their properties depend on the wind speed, duration and fetch length. Higher wind 
speeds, longer durations and longer fetches all result in larger wind-generated waves. 
However, in the nearshore environment, wave properties are also affected by interactions 
with bathymetry, tides and surface currents. These nearshore features can alter a wave 


field’s amplitude, speed and direction of propagation. 


The primary influence on nearshore wave behavior is the effect from bathymetry. 
As deep-water waves propagate into the nearshore environment, the water depth 
decreases and begins to affect the propagation speed, amplitude and direction. This effect 
is most evident near the beach where shoaling bathymetry slows down the approaching 
waves causing the wave amplitude and steepness to increase resulting in the breaking 
wave observed in the surf zone. The variable nearshore bathymetry also interacts with 
wave propagation by refracting the wave crest, altering the wave’s direction of 
propagation. This refraction process typically causes a focusing of wave energy on sand 


bars and coastal headlands. 


A second influence on nearshore wave behavior is that of tides. Ocean tides are 
most noticeable in the swash zone on the beach where the water level moves forward and 
backward throughout the day. Tides are described by three categories based on their 
frequency and variability: diurnal, semidiurnal, or mixed (see Figure 1). The frequency 
and characteristics of the tide are controlled by a combination of harmonic tidal 
constituents and vary based on location. Mixed tides, most common on the west coast of 
the U.S., experience two high tides and two low tides per day, with one of each being 
larger than the other. The large tidal movement of water into and out of an inlet can drive 
strong currents. Ebb currents (flowing out of a bay) and flood currents (flowing into a 


bay) can be extremely strong, particularly in narrow tidal inlets depending on the tidal 


range and the size of the bay. These tidally-generated currents can affect the propagation 
speed of a wave, thus altering the wave properties in a manner similar to the shoaling and 


refraction effects of variable depth. 
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Figure 1. Examples of tide characteristics (Fig (2.18a) in Open University 1999). 


In addition to the characteristics of the daily fluctuations, relative tidal strengths 
also observe a cyclical oscillation based on the lunar cycle. Spring tides occur twice per 
month coinciding with the celestial alignment of the sun, earth and moon. This results in 
the largest intra-daily differences in high tide and low tide. These large tidal differences 
result in the strongest tidally-generated currents. Conversely, neap tides occur when the 
sun, earth and moon are offset by 90° resulting in the smallest intra-daily difference 


between high tide and low tide, and accordingly weaker tidally-generated currents. 
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Other surface currents also affect nearshore wave behavior. Nearshore wind and 
wave-driven currents typically flow alongshore and generate eddy circulations and rip 
currents depending on local bathymetric features. These surface currents and circulations 


interact with approaching wave fields. 


The variability in direction and strength of currents can have direct impacts on the 
evolution and propagation of a wave field. An opposing surface current will tend to 
amplify wave heights at focal regions due to the refraction of wave energy about the 
current (Wolf and Prandle 1999). This situation can result in variable wave heights and 


directions in regions of strong wave-current interactions. 


These combined effects of bathymetry, tides and surface currents in the nearshore 
environment can result in large and abnormal wave conditions making ocean wave 
measurements difficult. Although traditional Eulerian measurements remain ideal for 
providing a robust time series of wave measurements at a specific location, they are not 
ideal for capturing the spatial variability which is often seen in the nearshore 
environment. A large array of fixed instruments or moored buoys would need to be 
deployed to capture such spatial variability, which is cost prohibitive. Additionally, heavy 
shipping channels pose a hazard to moored buoys where the buoys might sustain damage 
and accessibility for maintenance is restricted. The moorings themselves can also be 
negatively impacted in a dynamic nearshore environment, particularly in areas with 
rapidly changing bathymetry or extreme surface currents which can put undue strain on 
the moorings. On the other hand, drifting wave buoys that follow along with the surface 
current are able to capture both the temporal and spatial variability across this dynamic 
environment. Recent advances in inexpensive instrumentation have enabled the 
deployment of a large array of cost-effective buoys that are capable of observing the 
temporal and spatial evolution of the wave field across the nearshore environment. These 
large arrays of drifter buoys have the potential to provide accurate and reliable wave 
measurements in this dynamic environment which are essential for the improvement and 


validation of nearshore numerical models. 


Over the past five decades, various models have been developed to characterize 
and predict ocean waves such as WAM (Wave Action Model) and WAVEWATCH III. 
These models invariably assume Gaussian and homogeneous statistics and perform well 
in deep-water environments where the wave statistics are relatively weakly nonlinear and 
slowly varying. For nearshore applications, wave models such as SWAN (Simulating 
Waves Nearshore) or coupled wave-current-transport models such as DELFT-3D have 
evolved to accurately capture wave-current and wave-bathymetry interactions. However, 
accurate field observations that are required to initialize and validate these models are 


scarce. 


B. APPLICABILITY OF LINEAR WAVE THEORY 


Numerical wave prediction models and techniques for the analysis of ocean wave 
observations are often based on the principles of linear wave theory. This theory 
simplifies the solution to the wave equation by assuming that the wave amplitude is much 
smaller than both the wavelength and the water depth. The first-order solution for a 
sinusoidal wave obeys a linear dispersion relationship relating the wave frequency (@), 


the wavenumber (4) and the water depth (#1): 
@ = gk tanh(kH) (1) 


This linearized theory also gives approximate solutions for wave-induced fluid 
motions at and below the surface, and the associated pressure field. The accuracy of 
linear wave theory is well established. Thornton and Krapohl showed that measured wave 
orbital velocity fluctuations can be related to sea surface height through linear wave 
theory (1974). Herbers et al. (1992) demonstrated that the linear transfer functions 
between wave orbital velocities and bottom pressure (wave height) are accurate to within 


a few percent. 


The robustness of linear wave theory has led the way for a variety of in-situ 
instruments and remote sensing platforms than can be used for observing ocean waves 


with varying degrees of accuracy, resolution, and spatial/temporal coverage. 


C. OCEAN WAVE SENSING STRATEGIES 


Early sensors include simple tide staffs that can be deployed at any location to 
observe the instantaneous sea surface height, but this is a much localized and manual 
measurement. Radar altimeters deployed onboard aircraft and satellites are capable of 
providing wave height measurements based on the differential reflection of a radar pulse 
from the crest and trough of a wave (Stewart 2008). The satellite-based altimeter is ideal 
for observing a large area over a long period of time, but does not resolve the detailed 
wave motion. Coastal surface-based radars have been used for measuring ocean waves 
based on Bragg scattering from small-scale roughness elements on the sea surface. 
Similar to altimeters, these sensors are ideal for long-period observations, but are 
restricted to local, coastal regions, and their wave height measurements are frequency- 
dependent (McWilliams 2005). Additionally, the transfer functions between the radar 
back scatter and the surface wave heights are less well understood. Pressure sensors have 
been developed that are capable of observing the pressure differential of passing waves. 
These bottom-mounted sensors are extremely portable and ideal for deploying in a 
desired location; however, they are restricted to shallow water depths in order to avoid 
the attenuation of the “wave-induced pressure fluctuations” (Stewart 2008). Similar to 
pressure sensors, ocean current meters have been developed that are capable of observing 
wave heights through the changing water particle velocities associated with passing 
waves, but these sensors also only provide localized measurements (Thornton and 


Krapohl 1974). 


Still, the most common sensor for measuring wave heights has been the wave 
buoy, which physically observes the change in sea surface elevation. Wave buoys can 
either be moored in a specific location or allowed to drift with the current. The first 
generation of wave buoys relied on heave, pitch and roll measurements in order to 
determine the wave height and sea surface slope (H/P/R buoys). These measurements 
were made by accelerometers integrated into a motion reference unit (MRU). The vertical 
accelerations (heave measurements) were used to calculate the wave height, while the 
pitch and roll measurements were used to calculate the sea surface slope. Through the 


combination of the “vertical accelerations and the slope of the water surface...the 
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directional spectrum [was] estimated” (Long 1980). A modern example of an H/P/R buoy 


is the NOAA discus buoy (see Figure 2) (Krogstad 2003). 





Figure 2. A NOAA 3m H/P/R discus buoy (From 
http://www.ndbc.noaa.gov/images/buoys/3m. jpg). 


Further improvements to H/P/R buoys introduced a magnetic compass for heading 
as well as advancements in the accelerometers. Similar to the H/P/R buoys, these 
displacement buoys measured the wave height from vertical accelerations; however, the 
wave direction was obtained from a combination of the compass heading and horizontal 
accelerations. A modern example of these buoys is the Datawell Waverider buoy (see 


Figure 3) (Krogstad 2003). 








Figure 3. A 0.9m Datawell Waverider displacement buoy (From 
http://www.datawell.nl/inhoud.php?id=3). 


The most recent improvement to wave buoys was the addition of Global 
Positioning System (GPS) receivers (Krogstad 2003). Although GPS accuracy has varied, 
this step provided a direct measurement of both the wave height and wave direction from 
the vertical and horizontal displacements of the buoy. Continued improvements to the 
accuracy of GPS receivers have provided increasingly accurate measurements of buoy 
displacements, allowing for long distance tracking of the buoys. Improvements have also 
been made to the GPS sensing strategy reducing the size and cost of the GPS receivers 
(Herbers et al. 2012). The development of inexpensive GPS drifter buoys shows promise 


for deployments in large numbers to map the evolution of waves and surface currents. 


D. GPS BACKGROUND 


The evolution and accuracy of the GPS system is central to its employment in 
modern wave buoys used in this study. The GPS System is a passive, satellite-based 
positioning system originally designed for use by the Department of Defense (DoD) and 
the U.S. military, currently consisting of 36 active positioning satellites. These satellites 


transmit a coded signal that is received by a receiver on Earth. The travel time of the 
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coded signal is used to determine the range to the satellite — referred to as pseudo- 
ranging. The receiver uses pseudo-ranges from multiple satellites in order to triangulate 
its position on Earth (see Figure 4). Since satellite signals require a line of sight, they are 
only visible in a 180° hemispheric viewpoint. This underlying geometry provides for 


greater accuracy in the horizontal plane than in the vertical plane (El Rabbany 2006). 





Figure 1.5 The basic idea of GPS positioning. 





Figure 4. The principle of positioning through trilateration 
(Fig (1.5) in El-Rabbany 2006). 


The coded signal transmitted by the GPS satellite contains one of two binary 
codes. Low-cost commercial GPS receivers are typically single-frequency receivers that 
are capable of only receiving the coarse-acquisition code (C/A-code). The length and 
repetition rate of this code results in a real-time positioning resolution of 300m. More 
expensive survey-grade GPS receivers are typically dual-frequency receivers which can 
receive both the C/A-code and the precision code (P-code). The length and higher 
repetition rate of the P-code allows for real-time positioning accuracies within 30m. Until 
2000, the U.S. Government intentionally introduced some error into the GPS code known 
as selective availability (SA) in order to prevent high accuracy real-time positioning by 
commercial and non-military GPS receivers. This function has since been deactivated and 


the accuracy of the C/A-code has returned to 300m (El Rabbany 2006). 


The GPS signals are used to for real-time positioning by one of two techniques: 


code-phase or carrier-phase positioning. The code-phase positioning is the less accurate 
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of the two methods and involves the comparison of the receiver-generated code with the 
satellite-generated code, which results in pseudo-range accuracies that are dependent on 
the length of the code (approximately 300m or 30m). Assuming that the messages were 
initiated at the same time, the time delay required to synchronize the two messages 
corresponds to the travel time of the transmitted message, thus providing the pseudo- 
range to the satellite (see Figure 5). There is some error in this calculation due 
atmospheric propagation as well as inherent timing issues between the satellite and the 


receiver (El Rabbany 2006). 


Carrier-phase positioning involves the measurement of the total number of 
sinusoidal oscillations of the carrier frequency during the transmission (see Figure 5). 
This results in pseudo-range accuracies that are dependent on the wavelength of the 
carrier signal (19cm for the L1 frequency); however, the signal is purely sinusoidal and 


therefore the total number of oscillations cannot be readily measured. 
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Figure 5. GPS positioning techniques. The left panel shows the concept of code-phase 
positioning. The right panel shows the concept of carrier-phase positioning 
(Fig (2.4) and Fig (2.5) in El-Rabbany 2006). 


In order to improve real-time positioning accuracy of code-phase positioning, 
various differential correction techniques can be implemented in order to reduce some 
systematic errors. The most common technique is the use of a space-based augmentation 
system (SBAS). GPS receivers located at known geodetic positions across the country 


(and the world) calculate corrected pseudo-ranges to the satellites based on their known 


position. These corrections are transmitted to a satellite, which are then broadcast to other 
GPS receivers on the ground. This improves the real-time accuracy of standard GPS 


receivers (300 meters) to less than tens of meters (El Rabbany 2006). 


An additional improvement to modern GPS receivers is their measurement of 
velocity. Speed used to be measured by positional differences which contained tens of 
meters of uncertainty. Just as carrier-phase positioning began utilizing the characteristics 
of the carrier frequency, GPS receivers began measuring the Doppler shift of this 
frequency in order to determine the receiver’s velocity. Although this doesn’t provide 
real-time position, these measurements result in very accurate velocity measurements 


(within 10cm/sec) as opposed to absolute positional accuracies of tens of meters. 


The improvements in low-cost GPS receivers, such as the Locosys GT-31, which 
utilize SBAS DGPS corrections and measure velocities through the Doppler shift of the 
carrier frequency, result in accurate ocean wave and current measurements. The 
employment of these types of GPS receivers on drifter buoys provides the ability to 
measure wave-induced surface motions anywhere in the world, particularly the dynamic 


nearshore environment (see Herbers et al. 2012). 


E. ACCELEROMETER BACKGROUND 


Just as GPS has advanced over recent decades, accelerometers have also evolved, 
providing greater accuracy at ever reducing size and cost. Originally large mechanical 
devices, accelerometers provide measurements of accelerations in specific directions 
based on the displacement of a mass. However, the development of the MEMS (micro- 
electro-mechanical systems) sensor has been instrumental in making accelerometers 


smaller, more sensitive, and reliable. 


MEMS accelerometers use capacitance measurements between two surfaces to 
determine the acceleration of an object. A movable proof mass with extended arms is 
suspended on a spring between arrays of fixed arms (see Figure 6). As the accelerometer 
is moved, the suspended proof mass moves up or down in relation to the fixed array. This 
motion changes the relative distance between the arms of the suspended proof mass and 


the arms of the fixed array. Voltage sensors on the arms measure capacitance. The 
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variable capacitance measurements between the arms of the fixed array and the arms of 
the suspended proof mass are used to calculate the displacement of the suspended mass, 


and thus the acceleration (Andrejasic 2008). 
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Figure 6. MEMS accelerometer design (From Andrejasic 2008). 


The fundamental design of MEMS accelerometers measures both static (due to 
gravity) and dynamic (due to motion) accelerations, which allows these sensors to be 
oriented in any direction (Andrejasic 2008). Modern MEMS accelerometers have sensor 
plates aligned orthogonally in order to measure accelerations in all three axes (x, y, Z). 
Since the sensors measure static as well as dynamic acceleration, the vertically-aligned 
axis on a multi-axis accelerometer can be identified by the measured gravitational 
acceleration. The reduced size and improved sensitivity of these sensors make them ideal 


for integration into wave buoys. 


F. RESEARCH MOTIVATION 


The effects of the nearshore environment on ocean waves can be extreme and 


result in abnormal wave conditions which are hazardous to shipping and maritime 
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operations. The effects of bathymetry and tidal currents around inlets have been shown to 
focus wave energy in specific locations, setting up strong and variable surface currents. 
Bathymetric features can cause wave heights to increase and waves to break resulting in 
potentially hazardous conditions. Strong currents can also affect safety of navigation and 
generate strong eddy circulations. The effects of focused wave energy and strong surface 
currents can combine to result in shifting sand bars, sand waves, and beach erosion. 
Coastal zone management agencies continually strive to understand and mitigate these 
effects in order to preserve coastlines, coastal habitats and civilian coastal infrastructure. 
Nearshore waves and currents are constantly monitored for potentially hazardous 


conditions that affect swimmers, shipping and small boat operations. 


The U.S. Navy as well as civilian maritime agencies, continue to require safety of 
navigation within this dynamic nearshore environment in order to transit in and out of 
desired ports and harbors. This requires an accurate understanding and improved 
prediction capability of ocean waves. Although the accuracy of ocean wave models 
continues to improve, their improvement relies on the continued validation of their output 
and increasing accuracy of the initialization conditions. Recent research has shown the 
promise of using GPS wave buoys for ocean wave measurements. In this study we further 
refine an inexpensive GPS-tracked buoy (Herbers et al. 2012). Specifically, we add an 
inexpensive accelerometer to the GPS drifter buoy which provides more robust and 
accurate vertical surface height measurements, while maintaining a cost-effective 
Lagrangian drifter buoy. This buoy design is advantageous for deploying in a large array 
in order to accurately measure the spatial variability of wave fields in a dynamic 


nearshore environment. 


12 


Ul. FIELD EXPERIMENT AND DATA 


The field experiment was conducted at two locations during four cruises onboard 
the R/V Point Sur. The first part of the experiment involved the continued development 
and validation of a new wave resolving drifter (WRD) buoy. The second part of the 
experiment involved the deployment of a large array of these WRD buoys in and around 
the entrance to San Francisco Bay, where strong tidal currents and an ebb-tidal shoal 


interact with the approaching ocean swell creating strong wave-current interactions. 


A. BUOY VALIDATION 
1. Experiment Location 


The validation of the WRD buoy was performed approximately 30nm northwest 
of Monterey, CA (see Figure 7). The experiment was conducted from the R/V Point Sur 
based out of Moss Landing, CA. This validation site, well offshore of Monterey Bay in 
an average depth of 1700m, was selected to ensure that the surface waves were not 
affected by any bathymetric or coastal sheltering effects. The strong (24 mph) and steady 
wind out of the northwest (325°) allowed for the measurement of a fully developed wave 
field that was relatively homogenous in space and time. Additionally, the location was 
approximately 11nm from the National Data Buoy Center (NDBC) Station 46042, which 


provided coinciding environmental conditions. 
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Figure 7. Drift track of the validation deployment offshore of Monterey Bay, CA. The 
position of NOAA Buoy 46042 is noted (From Google Earth). 


The validation dataset was collected on April 29th, 2012 during one continuous 
eight-hour deployment. A total of four wave buoys were deployed: one Datawell 
Waverider buoy and three WRD buoys. The Datawell was deployed for comparison and 
treated as a “ground truth.” The buoys were deployed simultaneously at the same location 
from the working deck of the R/V Point Sur, ensuring that all four buoys observed the 
same wave field. During the eight-hour deployment, the Datawell and WRD buoys 
experienced slightly different drift characteristics, and were separated by 2km by the end 
of the deployment. 


2. Atmospheric Conditions 


Local atmospheric conditions were recorded by NDBC Station 46042 which is a 
3m discus buoy maintained by the National Data Buoy Center. This station provided 
hourly observations of atmospheric pressure, air temperature, surface wind speed (at 10m 
above the sea surface), surface wind direction, wave height, dominant wave period, 


average wave period, mean wave direction and sea surface temperature. 
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During the validation experiment, surface winds were out of the northwest (327°) 
at approximately 11m/s. The average wave height was 2.8m coming out of the northwest 
(310°) with a dominant wave period of 10 seconds. Hourly observations from Station 


46042 on April 29, 2012 are listed in Table 1. 


Table 1. | Atmospheric and wave conditions reported by NDBC Buoy 46042 during 
the period of the validation deployment. 


Pf are | Ma toed | i | te | 
(UTC) (m/s) (m/s) (m) (sec) 
| Apr29_ | oos0 | 327, | 94 | | 314 | 2.96 


| | 050 | 323 | 10.7 | 125 | 306 | 29 | 11.43 | 
| so | 320 | tas | 144 | 306 | 28 | 12.12 
| 10.8 


| 0350 | 323 | 13.2 | 306 | 278 | 11.43 
| | 450 | 319, | 109 =| 133 | 320 | 272 | 10 
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Figure 8. Wind speed and offshore wave height recorded at NDBC Station 46042 for 
the validation deployment. 


Ds Datawell DWR-G Buoy Sensors and Configuration 


The Datawell was deployed as the accepted industry standard for GPS wave 
buoys. The model that we employed was the 0.4m diameter Datawell Waverider — GPS 
(DWR-G4), which retails for approximately $14,000 (see Figure 9). This buoy is 
constructed from stainless steel, capable of withstanding a collision and corrosion, and is 
specifically designed for the marine environment. The buoy utilizes Doppler velocity 
measurements from an integrated GPS receiver to calculate horizontal and vertical 
excursions. The GPS receiver is capable of resolving positional displacements of 0.01m 
and direction within 1.5°. The GPS receiver operates at a sampling frequency of 1.28Hz. 
The onboard processor utilizes a proprietary algorithm, which converts the observed 
horizontal and vertical position displacements into wave height spectra and directional 


spectra. An onboard data logger records both horizontal and vertical position 
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displacements as well as the computed wave height and direction spectra. The DWR-G 
records one full GPS fix and transmits its position once every 30 minutes (Datawell 


2010). 


For this experiment, the DWR-G was also outfitted with the GT-31 GPS receiver 
and X6—2 accelerometer described in the next section for the purpose of sensor 


comparison and validation (see Figures 9 and 10). 





Figure 9. Datawell Waverider Buoy (DWR-G) drifting off of San Francisco Bar (left 
panel). DWR-G drifting off of Monterey Bay, CA (right panel). 





_ 


Figure 10. DWR-G external sensor configuration (left panel). GT-31 GPS receiver is 
visible on the bottom left. Two GoPro cameras are attached to the top of the 
buoy. DWR-G internal sensor configuration (right panel). The blue X6—2 
accelerometer is visible in the top left of the internal board, while a Yost 
accelerometer (for future testing) is visible in the center. 
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4. Wave Resolving Drifter (WRD) Buoy Sensors and Configuration 


The WRD buoy developed at NPS is built around a 0.4m diameter Jim Buoy float 
and costs less than $1,000 per buoy. The WRD is outfitted with a Locosys GT-31 
handheld GPS receiver and a Gulf Coast Data Systems X6—2 three-axis accelerometer. 
The buoys are also outfitted with a NaviSafe marine light, Merlin Systems radio 
transmitter and Garmin DC-40 GPS tracking collar for real-time position tracking (see 
Figure 11 and 12). The GT-31 GPS receiver is waterproof; however, it is secured inside a 
clear plastic bag for added protection in the nearshore environment. The GT-31 is secured 
near the top of the WRD in order to ensure optimal line-of-sight in large waves, 
minimizing satellite drop-outs. The X6—2 accelerometer is secured inside a watertight 
Otterbox case, which is attached approximately 0.3m beneath the buoy. The Otterbox 
case is suspended between the buoy and 20lbs of anchor chain used to ballast the buoy 
(see Figure 15). The anchor chain serves three purposes: it provides adequate ballast in 
order to ensure the centerline of the buoy rides at the water surface, it ensures the 
accelerometer remains consistently vertical while minimizing tilt, and it increases the 
submerged surface area of the buoy on which the current can act, to keep it well coupled 


to the surface drift. 





Figure 11. Wave Resolving Drifter (WRD) drifting inside Golden Gate Channel (left 
panel). WRD buoys drifting off of Monterey Bay, CA (right panel). 
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Figure 12. WRD Buoy. GT-31 GPS receiver (front), Garmin GPS tracking collar (left), 
Merlin RF transmitter (right) and NaviSafe marine light (back) are mounted 
externally on the buoy. 


a. Locosys GT-31 GPS Receiver 


The Locosys GT-31 is a single-frequency, handheld GPS receiver that 
costs approximately $135 (see Figure 13). This GPS receiver was previously used on 
earlier versions of GPS drifters developed for offshore wave observations (Herbers et al. 
2012). The unit utilizes a built-in patch antenna and lithium-ion polymer battery making 
it compact and portable. The battery is capable of 41 hours of continuous operation 
allowing it be used for day-long deployments. The receiver operates on the L1 carrier 
frequency, utilizing the C/A-code for pseudo-ranging. Additionally, the GPS receiver 
receives DGPS corrections via a space-based augmentation system (SBAS) allowing for 
greatly improved pseudo-ranging and positioning. The sampling frequency of the 


receiver is 1Hz (Locosys 2007). 
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The GPS receiver records measured horizontal and vertical displacements 
in meters with an uncorrected accuracy within ten meters. The receiver also measures the 
Doppler shift of the carrier frequency providing velocity measurements accurate within 
0.1m/s (Locosys 2007). Raw data is stored internally on a 2GB SanDisk SD card 


allowing for autonomous operation. 





Figure 13. GT-31 GPS receiver. 


b. X6-2 Accelerometer 


The Gulf Coast Data Concepts (GCDC) X6—2 is a_ three-axis 
accelerometer that costs approximately $115. The accelerometer is extremely small, 
slightly larger than an average thumb drive and weighs only 1.30z (see Figure 14). The 
accelerometer contains an internal lithium-polymer battery that is capable of 36 hours of 
continuous operation at the highest sampling frequency (320Hz) allowing it to be used 
for day-long deployments. The sampling frequency is adjustable and set to 10Hz for this 
experiment, providing much longer battery life. The accelerometer is capable of 
measuring accelerations up to 6g with an accuracy of 0.004g in the horizontal (X and Y 


directions) and 0.006g in the vertical (Z direction) (Gulf Coast Data Concepts 2010). 
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Figure 14. X6—2 accelerometer. 





Figure 15. Tethered configuration of the X6—2 accelerometer. The X6—2 is secured inside 
a waterproof Otterbox and is ballasted by 20lbs. of chain link. 
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c. Real-time Position Tracking 


Three different devices were implemented in order to aid in the real-time 
position tracking of the WRD buoys. Every WRD buoy was outfitted with a NaviSafe 
marine light, which cost approximately $70. The marine light was installed in order to aid 


in visibility and location of the buoys at night. 


Every WRD buoy was also outfitted with a small Merlin 1/3N MX 
transmitter that costs approximately $160 (see Figure 16). This transmitter was enclosed 
in a waterproof case that was sealed into the top of the WRD buoy. Each transmitter 
emits a unique radio signal between 216 and 219MHz at one second intervals. The MS- 
1000 Multi-Channel receiver was used for radio direction finding, providing coarse 


tracking and recovery. 
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Figure 16. The left panel shows the Merlin 1/3N MX radio transmitter. The right panel 
shows the MS-1000 Multi-Channel receiver used for radio direction finding 


(From http://www.merlin-systems.com). 


The Garmin DC-40 GPS tracking collar was attached to many WRD 
buoys and costs approximately $200. The GPS unit is relatively small, lightweight and 
water resistant (see Figure 17). The unit uses a built-in lithium-ion battery that is capable 
of continuous operation for up to 48 hours. This GPS tracking unit has a line-of-sight 
range of up to 11km (Garmin 2011). The Garmin Astro handheld device was used to 
monitor and track up to ten DC-40 GPS units per handheld device. 
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Figure 17. Garmin DC-40 GPS tracking collar. 


Table 2 outlines the sensor configuration, sampling frequency and 


observation accuracy for the two types of buoys that were deployed. 


Table 2. Wave buoys with sensor configurations and accuracies. 


Datawell DWR-G4 Wave Resolving Drifter (WRD) 
Datawell GT-31 X6-2 GT-31 X6-2 
GPS GPS Accelerometer GPS Accelerometer 


1.0 


Horizontal acceleration (m/s?) 
Vertical acceleration (m/s?) 
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B. NEARSHORE WAVE EVOLUTION 
i Experiment Location 


The nearshore buoy deployments were conducted during four cruises in and 
around the Golden Gate channel (see Figure 18), at the entrance to San Francisco Bay, 
CA. The San Francisco Bay is the largest coastal bay/estuary on the west coast of the 
United States and is characterized by large broad shallow expanses with a relatively 


narrow and deep navigable channel (Conomos et al. 1985). 


The Golden Gate Channel is situated at the mouth of San Francisco Bay. It is 
bounded by the San Francisco peninsula to the south and the Marin peninsula to the 
north. This channel is the only connection between San Francisco Bay and the Pacific 
Ocean. The width of the Golden Gate channel varies from 1.6 km to 4.0 km, with a 
maximum depth of 113m due to the hydrodynamic scouring from the intensified water 


flow through the constrictive geography. 
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Figure 18. Location of the nearshore deployment near San Francisco, CA. Example drift 
tracks following the ebb current are shown. The position of NDBC Buoy 
46026 is noted (From Google Earth). 
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The combined tidal freshwater river inflow to Suisun Bay in the northern reach of 
San Francisco Bay results in an “enormous tidal prism of 2 x 10° cubic meters (528 
billion gallons)” (Barnard et al. 2006). This tidal prism produces extreme tidal currents 
into and out of San Francisco Bay where peak flood currents (flowing into the bay) can 
reach speeds of 2.5 m/s and peak ebb currents (out of the bay) can reach speeds of 3.0 
m/s. This large tidal prism also results in extreme sediment transport. Sediments washed 
down from the Sierra Nevada Mountains through the Sacramento and San Joaquin Rivers 
are carried out of San Francisco Bay through the Golden Gate Channel. Outside the 
channel these sediments have been deposited over time resulting in a crescent-shaped 
ebb-tidal shoal that is only five meters deep at its shallowest location. This ebb-tidal 
shoal surrounds the mouth of the Golden Gate Channel and presents a significant and 


unique bathymetric feature that affects the approaching wave field (see Figure 19). 
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Figure 19. Ebb-tidal shoal located at the mouth of the San Francisco Bay 
(From Elias and Hansen 2012). 


In addition to the ebb-tidal shoal, the extreme tidal currents and sediment 
transport through the Golden Gate Channel result in large sand waves at the mouth of the 


channel (see Figure 20). These sand waves have an average wavelength of 82 meters and 
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an average height of six meters (Barnard et al. 2006). The sand wave field at the mouth of 
the Bay is comprised of more than 40 individual sand waves and occurs in depths as 
shallow as 30 meters. Although the largest of these sand waves occur in the deeper 
portions of the channel and do not pose any hazard to shipping, they do influence the 


water depth, ultimately affecting the surface current and wave evolution. 
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Figure 20. Sand waves located near the mouth 


2. Atmospheric Conditions 


NDBC Station 46026 is a 3m discus buoy maintained by the NDBC. This station 
provides hourly observations of atmospheric pressure, air temperature, surface wind 
speed (at 10m above the sea surface), surface wind direction, wave height, dominant 


wave period, average wave period, mean wave direction and sea surface temperature. 


During the February deployment, offshore winds came out of the northwest 
(317°), increasing overnight on February 15th before abating late in the morning on 
February 16th. The surface winds peaked at 13.8m/s at 1950. The offshore wave field 
steadily came out of the west-northwest and varied in height between 2.14m and 3.26m. 
Hourly observations from station 46026 for February 15-16 are given in Table 3 and 
Figure 21. A peak ebb current of 4.2kts occurred in the Golden Gate channel at 
1718UTC. 
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Table 3. Atmospheric and wave conditions reported by NDBC Buoy 46026 for 
February 15th-16th, 2012. 


Time Wind Spd} Gust Wave Het | Wave Per| Predicted 
(UTC) (m/s) (m/s) (m) (sec) Current 
Tk 


| 71 | 88 | 284 | 275 | 16 | 
ee 
a ee eae ee ee 
| 1650 | 312 |e | 154 | 278 | 246 | 14.81 | Ebb 
| aso | 3aa_ | 13.3 | 162 | 282 | 252 | 13.79 | Ebb(4.2) 
| | aso | 312 | 13.3 | 167 | 286 | 266 | 1481 | Ebb 
SEITEN ae aEnenEe Ueeeney Camere eee eoaela a 
| 320 | 121 | 14.4 | 288 
ee 
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Figure 21. Wind speed and offshore wave height recorded at NDBC Station 46026 for 
the February deployment. 


During the April deployment, offshore winds were relatively consistent coming 
out of the northwest at 10m/s. The offshore wave field steadily came out of the west- 
northwest and slowly increased in height during the deployment to a significant wave 
height of 2.56m. Hourly observations from station 46026 on April 27 are given in Table 4 
and Figure 22. A peak ebb current of 3.5kts occurred in the Golden Gate channel at 
1S558UTC. 
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Table 4. | Atmospheric and wave conditions reported by NDBC Buoy 46026 for April 
27th, 2012. 


Time Wind Spd | Gust Wave Hgt | Wave Per | Predicted 
(UTC) (m/s) (m/s) (m) (sec) Current 
(kts) 


|Apr27_ | 0050 | 292 | 92 | 107 | 268 | 106 | 16 | Flood | 
| nso | 295 | 78 | 93 | 247 | 135 | 1481 | Slack | 
| torso | 293 | 2 | 7.3 | 297 | 132 | 10 | Ebb 
| | aso | 296 | 74 | 88 | 294 | 137 | 9.09 | Ebb(16) | 
| | 450 | 295 | 8 | 77 | 303 | 15 | 909 | Ebb 
| sf 550 | 298 | st | 7.4 | 288 | 165 | 909 | Ebb 
| | 650 | 300 || S| 76 | 284 | 166 | 1081 | Slack | 
| if o7so | 299 | 5.4 | 67 | 290 | 163 | 833 | Flood | 
| | aso | 296 =| | ss | a9 | 164 | 10 | Flood 
| o9s0 | 298 | 69 | 79 | 283 | 166 | 10 | Flood (15) | 
| | 1050 | 295 | 79 | 95 | 293 | 158 | 10 | Flood | 
| Cf tt50 | 299 | 8s | 95 | 294 | 154 | 10 | Flood 
| 250 | 305 | os | a3 | 249 | 71 | 13.79 | Slack | 
| 350 | 313, | 8s | 101 | 299 | 195 | 10 | Ebb | 
| Ct 1450 | 310 | So | tos | 295 | 183 | 909 | Ebb | 
| sso | ais | as | ans | 300 | 209 | 10 | bb(3.5) | 
| | eso | 316 | 99 | 119 | 302 | 207 | 10 | ebb 
| dt 750 | 316 | 106 =| 123 | 304 | 216 | 10 | ebb | 
| sf 1850 | 310 | = 108 | 129 | 291 | 219 | 10 | Ebb | 
| it i950 | aaa | aaa | 128 | 297 | 22 | 1081 | Slack | 
| | 2050 | 31a | 10.3, {| 123 | 300 | 228 | 1081 | Flood | 
| CT 250 | 3a7_| 104 | 123 | 297 _ | 256 | 1212 | Flood | 
| | 2250 | 320 | 0.4 | 122 | 279 | 2.27 | 13.79 | Flood (2.5) | 
| aso | 3a7_ | 105 | 124 | 298 | 24 | 1212 | Flood | 
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Figure 22. Wind speed and offshore wave height recorded at NDBC Station 46026 for 
the April deployment. 


3. Buoy Deployment 


On February 15th, 2012, 20 WRD buoys were deployed from R/V Questuary 
underneath the Golden Gate Bridge. The buoys were deployed individually at 
approximately 12 minute intervals. The buoys drifted out the channel on the outgoing ebb 
tide where they were later recovered at various times by R/V Point Sur. Some of the 
buoys remained adrift for more than 24 hours. Seven of the buoys were ultimately lost 
and/or damaged after washing ashore north and south of the channel. This deployment 


captured the broad circulation characteristics of the tidal inlet. 
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On April 27th, 30 WRD buoys were deployed from R/V Questuary in the center 
of the ebb-tidal jet (as determined by ADCP transects), underneath the Golden Gate 
Bridge. The buoys were deployed in groups of five. There was a one minute interval 
between buoys and five minute intervals between groups. The buoys were deployed on an 
increasing ebb-tidal current and drifted out the Golden Gate. The buoys were recovered 
by R/V Point Sur after they had drifted seaward of the San Francisco bar (ebb-tidal 
shoal). 
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Hil, DATA PROCESSING AND ANALYSIS 


A. DATA CONSOLIDATION 
1. Data Extraction 


The raw data was extracted from the DWR-G buoy and WRD sensors in order to 
create one consolidated data structure for analysis and comparison (see Table 5). The 
DWR-G buoy recorded raw horizontal and vertical displacements in a proprietary 
Datawell binary format. In addition, it recorded the internally processed wave height and 
directional spectra that were not used in this study. The raw data extracted from the 
DWR-G provided four data strings: GPS time, easting displacement, northing 
displacement and vertical displacement. These data strings were all sampled at 1.28Hz. A 
30 second moving average was removed from the three displacement data strings in order 


to remove any low-frequency signal (Herbers et al. 2012). 


The GT-31 GPS receiver recorded the GPS time, latitude, longitude, positional 
displacements and horizontal Doppler velocities on an internal SD card. The data 
extracted from the GT-31 provided six data strings: GPS time, easting displacement, 
northing displacement, vertical displacement, easting Doppler velocity and northing 
Doppler velocity. All of these data strings were sampled at 1Hz. Similar to the DWR-G 
time series, a 30 second moving average was removed from the displacement data strings 
in order to remove any low-frequency trend or signal (Herbers et al. 2012). All of the 
DWR-G and GT-31 data stings were linearly interpolated to 10Hz in order to match the 


sampling frequency of the accelerometer. 


The first step of the consolidation aligned all of the raw data strings on a 
synchronized time base. Since the DWR-G displacements were achieved through an 
integration of the observed GPS signal Doppler shift, the DWR-G time was internally 
interpolated between the buoy’s 30-minute GPS fixes (Datawell 2010). Therefore, the 


more accurate GT-31 GPS time was used as the absolute time. A cross correlation 


a3 


analysis between the easting displacements of the GT-31 and the DWR-G provided the 
timing adjustment required to align the DWR-G displacements in accordance with the 


GT-31 time (Herbers et al. 2012). 


The data extracted from the X6—2 accelerometer provided three data records: 
acceleration in the x-axis direction, acceleration in the y-axis direction and acceleration in 
the z-axis direction). The vertically oriented axis of the accelerometer was easily 
identified based on the large gravitational signal. The vertical acceleration string was 
integrated once in order to calculate the vertical velocity, and similar to the GPS time 
series, a 30 second moving average was removed in order to remove any low-frequency 
noise. These velocity data were then linearly interpolated to 10Hz. Next, a cross 
correlation analysis between the vertical velocities of the X6—2 and the easting and 
northing displacements of the GT-31 provided the timing adjustment required to align the 
X6—2 velocities in accordance with the GT-31 time. The vertical velocity string was then 
integrated once more (again removing a 30 second moving average) in order to calculate 
the vertical displacement. A final cross-correlation analysis between the vertical 
displacements of the X6—2 and the vertical displacements of the DWR-G verified that the 


X6-—2 displacements where aligned in accordance with the GT-31 time. 


Finally, the synchronized data strings from the sensors were consolidated into 
one data structure which contained: GPS time, GPS latitude, GPS longitude, DWR-G 
easting displacement, DWR-G northing displacement, DWR-G vertical displacement, 
GT-31 easting displacement, GT-31 northing displacement, GT-31 vertical displacement, 
GT-31 easting velocity, GT-31 northing velocity, X6—2 vertical acceleration, X6—2 
vertical velocity and X6—2 vertical displacement. All data strings were now sampled at 
10Hz and aligned on the GT-31 GPS time (see Table 5). This process allowed for the 


direct comparison of the time series from the various sensors. 
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Table 5. | Raw and consolidated data output from the various buoy sensors. 


ae 
Sample rate | Time source 
GT-31 time (UTC) GT-31 
GT-31 latitude GT-31 
GT-31 longitude GT-31 
Datawell time (UTC) [4.28 | Datawell | 100 | GT-31 
Datawell easting displacement GT-31 
Datawell northing displacement | 1.28 | Datawell | 10.0 | GT-31 
Datawell verticaldisplacement | 1.28 | Datawell | 10.0 |  GT-31 
0 GT-31 
| 10) | GT31_—| 10.0 31 
GT-31 vertical displacement | 10 ©=6| ~ ~=6GT-31—C | 10.0" GT-31 
GT-31easting Dopplervelocity | 1.0 | GT-31 | 100 ~ | GT-31 
GT-31 northing Dopplervelocity | 1.0 | GT-31 | 100 ~ | GT-31 
X6-2 easting acceleration | 100 | x62 | 100 | GT-31 
[X6-2northing acceleration _—|__—10.0_—s|_—x6-2_—|— 10.0 | GT-31_— 
[X6-2verticalacceleration _—|__—10.0_~—s | x6-2_— | 10.0 | GT-31__— 
[X6-2verticalvelocity | ~~ | | 10.0] GT-31_ 
[X6-2 vertical displacement |] | | 10.0 | GT-31__ 


GT- 


.0 
.0 
GT-31 northing displacement 


1 
1 
1 
GT-31 easting displacement 
1 
1 
1 
a 





2. Time Series Characteristics 


The consolidated time series of the validation dataset included 307,199 data 
points, corresponding to 8.53 hours of data. The nearshore datasets varied in length from 


1.25 hours to 25.38 hours of data depending on the buoy deployment. The Nyquist 


1 : 
frequency of the buoy’s sensors dn ==) determines the maximum frequency that 


can be used in the spectral analysis. The Nyquist frequency of the Datawell buoy is 
0.64Hz, while the Nyquist frequency of the GT-31 GPS receiver and the X6-—2 
accelerometer are 0.50Hz and 5Hz, respectively. So, the DWR-G is capable of resolving 
waves with a period as short as 1.56 seconds, while the WRD buoy is capable of 
resolving waves with a period as short as 2 seconds. The interpolation of the data strings 
to 10Hz was performed on the collected data, thus it does not affect the sensor’s Nyquist 


frequency, and does not increase the buoy’s ability to resolve shorter period waves. 
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3. Linear Wave Theory Transfer Functions 


The present dataset provides robust measurements of various properties of the 
wave orbital motion. Through the application of linear wave theory, these measurements 
can be used to estimate wave frequency and direction spectra, as well as bulk wave 


statistics (Herbers et al. 2012). 


Assuming that the wave height is sufficiently smaller than the wavelength and 
that the current speed (buoy drift rate) is sufficiently smaller than the wave propagation 
speed, linear wave theory relates the orbital motion of the water particles to the sea 


surface elevation (wave height) (Herbers et al. 2012). To derive these transfer functions, 


the surface elevation (7) is expressed as a sum of Fourier components with different 


frequencies (a) and propagation directions 





1( x,t) = Ay exp|i{k-¥-ar) ]+ ce. (2) 


where k is the wave-number vector and c.c. is the complex conjugate. The wave number 


magnitude k = | is related to w through the linear dispersion relation (Equation 1). 


Transfer functions are first derived in order to relate the vertical motion (velocity 
and acceleration) to the surface height (7) at a fixed location (x = 0) , neglecting the drift 
of the buoy. The vertical velocity w(t) follows from taking the first time-derivative of 7 


at x =0. Similarly, the vertical acceleration a(t) equals the second derivative of the sea 


surface displacement equation. 


w(t)= >) > -i@A,, exp(—iat)+c.c. (3) 
o 0 


a(t)=)°)-o’4,, exp(-it) +c.c. (4) 
o 0 


Next, transfer functions are derived in order to relate the horizontal motion 
(velocity and displacement) to the surface elevations (vertical displacement). The 


momentum equations relate the horizontal accelerations to the slopes of the sea surface: 
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du on 


dt” x oe 
ce! 
ao 6) 


Differentiating the sea surface displacement with respect to x and integrating 


with respect to t (atx =0) yields the horizontal velocities (u, v): 


u(t)= ps cos 04, exp(—i@t)+c.c. (7) 
v(t) “Dye sin 0A, exp(—iat)+c.c. (8) 


Finally, the horizontal displacements (x, y) are derived from integrating the 


horizontal velocities (Equations 7 and 8) with respect to ¢ once more, 


x(t) = Ty cos 0A, exp(—iat)+c.c. (9) 
o 9 O 
igk . 
y(t) =D D7 sin Ag exp(—iat)+c.c. (10) 
o @ 


The frequency-dependent coefficients of the transfer functions are summarized in 


Table 6. 


Table 6. Derived coefficients for the transfer functions used to relate measured 
surface motion to sea surface displacement. 


Vertical displacement 
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B. SPECTRAL ANALYSIS TECHNIQUES 


The time series of wave-induced surface motions can be analyzed spectrally 
through the use of the Fast Fourier Transform (FFT), which relates the time-dependent 
motions to the wave energy spectrum (£(f)). This analysis provides an assessment of how 
the wave energy is distributed across a specific frequency band. For example, a spectral 
peak can identify a dominant wave period of a dataset. A cross-spectral analysis was 


performed using various combinations of the horizontal and vertical motion time series to 
estimate the mean wave direction, 0 ( - ) and the directional spread, o, ( f ) as a function 


of frequency. Finally, the bulk wave statistics (significant wave height, mean wave 
period, mean wave direction and directional spreading) were estimated from the 


frequency and direction spectra. 


Bulk values of @ and o, are estimated by replacing the auto- and cross-spectra in 


Table 7 with the integrals over the same frequency range f; — fo. The auto-spectra, E,(f), 
E\,(), E-(f), and cross-spectra, Cy(f), Cy-(f), C\-(), for each surface motion time series 
were calculated utilizing the Matlab cross power spectral density function (CPSD). The 
CPSD function utilizes a specified FFT record length with a 50% overlap of samples. The 
resultant spectrum was then merged across a number of frequency bands. Next, following 
Herbers et al. (2012), the frequency was restricted to 0.05Hz — 0.3Hz in order to only 


examine the swell-sea frequency range (wave periods 3.3 — 20 seconds). 


1. Directional Estimates 


The directional Fourier coefficients (see Long, 1980 for definitions) were 
estimated from the auto- and cross-spectra using the relations in Table 7 which 
incorporate the transfer functions of Table 6 for the various measurement combinations. 
The first-order coefficients a;, b; depend on both vertical and horizontal wave-induced 
surface motions; the second-order coefficients a2, by depend only on horizontal motions 
(and thus are not affected by errors in the acceleration measurements, see e.g. Herbers et 


al. 2012). 
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Table 7. Estimates of the directional Fourier coefficients for different combinations 
of measurements. 





a1 





bi 


a2 





b2 














The mean wave direction (@ ) and directional spread (o,) were estimated from 


both the first- and second-order Fourier moments (a;, b;) (Herbers et al. 2012, Kuik 
1988). The first-order estimates, given by 


tan(9) = 
a 
: (11) 
O,= 2(1- Jai +b 
(12) 
are widely used in routine analysis of wave buoy observations. 
The second-order estimates 
tan (20 ) = 
7 (13) 
(1 ae ee 
Oo, =4\/——_—— 
: (14) 


suffer from a 180 ambiguity, but do not rely on vertical surface motions, and thus provide 


a useful consistency check on the quality of the directional wave data (Herbers et al. 
2012). 
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2. Wave Parameters 


The significant wave height and wave period are readily estimated from the wave 
energy spectrum (Herbers et al. 2012). The significant wave height is commonly defined 
as the average of the largest one third of the waves. Spectrally, this is equivalent to four 


times the variance of the wave energy spectrum. 


ht 
H,=4 || E(f)af (15) 
A 


where f; and fo specify the frequency range of interest. 
The mean wave period was estimated spectrally from the frequency-weighted 


mean of the wave energy spectrum. 


[e(na 
f= a (16) 
| f-E(N)af 
fi 


J: Wave Motion Analysis 


As noted before, wave frequency and direction spectra (and associated bulk 
statistics) were calculated using various combinations of time series of horizontal and 
vertical motions to compare the calculated statistics from different buoy sensors. The first 
processed method utilized both the horizontal and vertical DGPS displacements (x,y,z). 
This method was used for analyzing the DWR-G and some of the GT-31 data. The 
second processing method utilized horizontal DGPS displacements from the GT-31 (x,y) 
and vertical accelerations from the X6—2 (a). Finally, the third processing method utilized 
horizontal Doppler velocities from the GT-31 (u,v) and vertical accelerations from the 


X6-2 (a). 
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IV. WAVE RESOLVING DRIFTER VALIDATION 


The validation of the WRD buoy focused on evaluating the sensor accuracy and 
buoy response in energetic open ocean wave conditions using the well tested Datawell 
DWR-G buoy as a “ground truth.” The first step of the validation was a comparison of 
measurements from WRD sensors (GT-31 GPS receiver and X6—2 accelerometer) 
collocated on a DWR-G buoy with the DWR-G measurements. The second step of the 
validation was an evaluation of the overall performance of the buoy design by comparing 
the measurements from the three WRD buoys with measurements from the DWR-G 


buoy, all drifting in close proximity to each other. 


A. SPECTRAL PROCESSING 


Prior to conducting the comparisons, the time series were analyzed in order to 
determine the optimal processing parameters. An analysis of various FFT record lengths 
and merging of frequency bands was done. The results of these analyses were used to 


process all of the spectra and statistics. 


1. FFT Record Length 


Various FFT record lengths were analyzed to determine the optimal frequency 
resolution while still resolving the temporal variations in wave conditions. These record 
lengths included 1024 data points (1.7 minutes of data), 2048 data points (3.4 minutes of 
data), 4096 data points (6.8 minutes of data), 8192 data points (13.7 minutes of data) and 
16384 data points (27.7 minutes of data). 


A comparison of the statistics calculated from the DWR-G displacement data (see 
Table 8) showed that the bulk wave statistics calculated from FFT record lengths of 
16384 samples and 8192 samples agreed well with the statistics calculated from the entire 
time series. All of the bulk wave statistics (calculated from both horizontal and vertical 
motions) differed by less than 1% and less than 2.5% respectively. Specifically, the 
significant wave heights differed by less than 0.02m for an FFT record length of 16384 
samples and differed by less than 0.09m for a record length of 8192 samples. The mean 
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wave direction and directional spread differed by less than 1° for both of these record 


lengths. 


The statistics calculated from the other FFT record lengths (4096, 2048 and 1024) 
also agreed well, differing by less than 8.5% when compared to the entire time series. Of 
note, the significant wave heights differed by less than 6.5% (0.21m) for an FFT record 
length of 4096 samples. These statistics improved slightly for an FFT record length of 
2048 samples, differing by less than 3.5% (0.09m). The mean wave direction differed by 
less than 1° for these record lengths while the directional spread differed by less than 


2.5°. The mean wave period was within 4.5% (0.37sec) for these record lengths. 


Table 8. | Bulk wave statistics for various FFT record lengths compared to the 
statistics calculated from the entire time series using the Datawell DGPS 


displacement data. 


Datawell GPS Displacements 


Vertical Orbital Motions FFT:1024 | iff | FFT:2048 | Diff 


Significant wave height (m) 


3.31 


3.38 


2.37% 3.38 2.21% 


FFT: 4096 


Diff 


FFT: 8192 


Diff 


FFT: 16384 





Mean wave period (sec) 


8.17 


8.28 


1.33 % 8.25 1.05 % 





Mean wave direction (deg) 


309.94 


310.33 


0.39° 310.36 0.42° 





Directional spread (deg) 


26.99 


26.90 


~-0.09° 24.76 -2.23° 


Horizontal Orbital Motions FFT:1024 | Diff | FFT:2048 | Diff 


Significant wave height (m) 


3.38 


3.58 


FFT: 4096 


Diff 
1.54 % 


FFT: 16384 
3.39 





Mean wave period (sec) 


8.32 


8.66 


1.09 % 


8.34 





Mean wave direction (deg) 


311.48 


311.71 


-1.10° 


310.97 





Directional spread (deg) 


20.65 


20.16 














0.01° 





20.77 





Scatter plots comparing the statistics of the DWR-G displacements and various 
FFT record lengths are shown in Figure 22 with statistics from vertical displacements 
plotted as squares and statistics from horizontal displacements plotted as circles. The 
plots show excellent agreement for all of the FFT record lengths except for a record 
length of 1024 samples. The best agreement is seen using record lengths of 16834 and 
8192 samples. 
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Comparison of FFT Record Lengths 
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Figure 23. Scatter plot comparison of the bulk wave statistics using various FFT record 
lengths of DWR-G DGPS displacements compared to statistics using the 
entire time series. Statistics calculated from vertical displacements are plotted 
as circles, while statistics calculated from horizontal displacements are 
plotted as squares. The top left panel is significant wave height, the top right 
panel is mean wave period, the bottom left panel is mean wave direction and 
the bottom right panel is directional spread. The dashed lines represent +/- 
5% from the statistic of the entire time series. 


A comparison of estimates based on GT-31 DGPS displacements versus GT-31 
Doppler velocities also showed excellent agreement for record lengths of 16384 and 8192 
samples with differences within 1.5% and 2.5% respectively (see Table 9). Statistics 
calculated from shorter record lengths of 4096 and 2048 were somewhat degraded (errors 


within 10%) with larger errors for a record length of 1024 samples. 
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Table 9. — Bulk wave statistics calculated for various FFT record lengths compared to 
the statistics calculated from the entire time series using various 
configurations of the GT-31 DGPS displacements, GT-31 Doppler velocity 
and X6-2 acceleration data. 


Datawell GT-31 DGPS Displacements 


|__Vertical Orbital Motions___|_Total_| FFT: 8192 | _biff_| 
a ro er oe 
Mean wave period (sec) 9.1 0.01% 

| Meanwavedirection(deg) ——|__ 309308 -0.95" | 
| Directional spread (deg) Tat Tat 0. 


Horizontal Orbital Motions FFT: 8192 | pitt | 
Significant wave height (m) | 3.39 2.35% 


Mean wave period (sec) 8.8 8.9 1.97% 
Mean wave direction (deg) 311 310 -0.34° 


Directional spread (deg) p22 Tote | 








Datawell GT-31 DGPS Displacements and X6-2 Accelerations 


|_Frt:8192 | oift_| 
Mean wave period (sec) 8.0 8.1 0.60 % 

| 309 | 309 | -0.53° | 
Directional spread (deg) | 34 | 35] taae | 





Horizontal Orbital Motions FFT: 8192 | Diff | 
Significant wave height (m) | 3.31 7 3.39 | 2.35% 


Mean wave period (sec) 
Mean wave direction (deg) 
Directional spread (deg) 





: 8.9 1.97% 
311 310 -0.34° 


ee 


ne oo 
Ld Go 


Datawell GT-31 Doppler Velocities and X6-2 Accelerations 
Vertical Orbital Motions Total | FFT: 8192 | Diff | 
Significant wave height (m) 3.50 | 3.56 | 157% | 
Mean wave period (sec) | 81 | 060% | 


Mean wave direction (deg) 
Directional spread (deg) 34 0.82° 

Horizontal Orbital Displacements Total | FFT:8192 | Diff | 
Significant wave height (m) 3.50 
Mean wave period (sec) a 8.7 1.88% 

Mean wave direction (deg) | 310 ~=|_ ona | 
Directional spread (deg) | 22 =| ~0.06" | 








ly 





PE 


In addition to the bulk wave statistics, the wave frequency spectra were calculated 
estimated from the DWR-G displacement data were compared for different choices of the 
FFT record length (see Figure 24). As expected, these results show the trade-off between 
degraded resolution for short FFT record lengths versus statistical scatter for long FFT 
lengths. FFT record lengths of 16384 and 8192 samples appear to resolve similar spectral 
features across the entire frequency range, as is the case for an FFT record length of 4096 
samples for most frequencies (greater than 0.08Hz). Spectra from record lengths of 2048 


and 1024 samples show degraded frequency resolution at frequencies below 0.2Hz. 


44 


Analysis of FFT Window Length 
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Figure 24. Wave frequency spectra calculated from various FFT record lengths of the 
DWR-G time series. The top panel uses vertical displacements. The bottom 
panel uses horizontal displacements. 


2. Frequency Band Merging 


Similar to the FFT record length analysis, various numbers of frequency bands (3, 
5, 7 and 9) were merged together during the calculation of the wave frequency spectrum. 
This was done in order to smooth out the spectrum and identify spectral peaks of the 
swell-sea and wind-sea bands. However, as discussed in Section 3.A.2, merging 
frequency bands together decreases the frequency resolution. Figure 25 shows the wave 
frequency spectra calculated using DWR-G displacements and an FFT record length of 
8192 samples. The minimum number of frequency bands merged together was three. This 
spectrum has the highest frequency resolution, but is difficult to accurately identify 
spectral peaks due to the increased number and variability of data points. Merging more 
frequency bands together resulted in a smoother spectral curve making the identification 
of spectral peaks more obvious. The optimal results are observed through the merging of 
seven frequency bands. This provides the best frequency resolution with the clear 


identification of the spectral peaks. 
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Analysis of Frequency Band Merging 
Datawell GPS 























Figure 25. Wave energy spectra calculated from the DWR-G time series using an FFT 
record length of 8192 samples and merging various number of frequency 
bands. The top panel uses horizontal displacements. The bottom panel uses 
vertical displacements. 


B. SENSOR COMPARISONS 
i. Linear Interpolation of Raw Data 


The various sensors all had different sampling frequencies, which were linearly 
interpolated to 10Hz to create the consolidated dataset. Statistical and spectral analysis 
was performed comparing the interpolation of the DWR-G data from 1.28Hz to 10Hz and 
the interpolation of the GT-31 GPS time series from 1.0Hz to 10.0Hz. 


Bulk wave statistics were calculated for each time series based on both vertical 
and horizontal displacements. These statistics are listed in Table 10. The interpolation of 
the DWR-G time series introduced less than 2% error into the calculated statistics. The 
significant wave height differed by less than 1% (0.02m) from the original statistics. The 
mean wave direction and directional spread also differed by less than 1% (1.2°) from the 


original statistics. The largest error introduced by the interpolation was observed in the 
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calculation of the mean wave period; however, it was still within 2% (0.2sec) of the 
original statistics. There was not any significant difference between the vertical and 


horizontal displacements. 


The interpolation of the GT-31 time series introduced slightly larger errors, but 
were still within 5% of the original statistics overall. The significant wave height differed 
by less than 2% (0.14m) from the original statistics. The mean wave direction and 
directional spread differed by less than 3% (approximately 1.5°) from the original 
statistics. Similar to the DWR-G, the largest error introduced by the interpolation was 
observed in the mean wave period, which differed by 5% (0.4sec) from the original 
statistics. Overall, there was slightly better agreement from the statistics calculated from 
the horizontal displacements than those calculated by the vertical displacements, which 


was likely caused by the decreased vertical resolution of the GT-31 GPS receiver. 


Table 10. | Bulk wave statistics calculated before and after the linear interpolation of 
the Datawell DGPS and GT-31 time series. 


Sampling frequency (Hz) 1.26 10.0 Diff 1.00 10.0 Diff 
Significant wave height (m) -0.46 % -1.70% 
Mean wave period (sec) 8.1 8.3 1.93 % 8.6 9.1 5.08 % 
Mean wave direction (deg) 310 309 -1.26° 309 308 -1.42° 
Directional spread (deg) 27 27 -0.17° 39 41 1.16° 
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Similar to the statistical comparison, a spectral analysis was performed comparing 
the wave frequency spectrum from the original and interpolated time series of both 
sensors. The spectral comparison of the DWR-G time series is shown in the top half of 
Figure 26. The interpolated spectrum agrees extremely well with the original spectrum 
identifying similar spectral peaks and spectral energy levels. The spectral comparison of 


the GT-31 time series is shown in the bottom half of Figure 25. Again, the interpolated 
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spectrum agrees very well with the original spectrum. The two spectra identify similar 
spectral peaks and spectral energy for frequencies less than 0.17Hz. For frequencies 
greater than 0.17Hz, the interpolated spectrum shows similar spectral characteristics; 
however, the spectral energy levels are slightly decreased due to the dilution of high 
frequency energy caused by the energy leakage beyond the original Nyquist frequency 


introduced by the linear interpolation (Pearman et al. 2013). 
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Figure 26. Spectral comparison of the linear interpolation of the time series. The left 
panels are the DWR-G wave spectra from vertical and horizontal 
displacements. The right panels are the GT-31 wave spectra from vertical and 
horizontal displacements. 


Overall, this analysis shows that the linear interpolation of the original time series 
does not significantly affect the bulk wave statistics or the wave energy spectrum, 
particularly in the swell-sea frequency band. The linear interpolation of the GT-31 seems 
to slightly degrade wave energy at higher frequencies, corresponding to waves with 
periods shorter than five seconds which is approaching one half of the sensor’s actual 


Nyquist frequency. 
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2. Horizontal Wave Motion 


The GT-31 GPS receiver provided two different techniques for measuring the 
horizontal surface motions: traditional DGPS positioning and Doppler velocity 
measurements. The quality of these two data sources was evaluated through comparisons 
of the time series, bulk wave statistics and wave energy spectra with the DWR-G 


estimates. 


An example of the horizontal easting displacement time series (see Figure 27) 
illustrates the excellent agreement between the GT-31 DGPS displacements and the 
DWR-G displacements (see also Pearman et al. 2013). The peaks of the signals track 


well. 
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Figure 27. Displacement time series of DWR-G and GT-31 easting displacements. 


The bulk wave statistics from the various horizontal surface motions are listed in 
Table 11 for comparison. Statistics from the GT-31 DGPS displacements agreed 
extremely well with the DWR-G displacements. The significant wave height differed by 
less than 2% (0.04m) and the mean wave direction was identical to that calculated by the 
DWR-G buoy. The mean wave period differed by approximately 6% (0.13sec) from the 
DWR-G buoy. The largest difference was observed in the directional spread, which 
differed by approximately 7% (~1.5°) from the DWR-G buoy. 
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Overall, the bulk wave statistics from the GT-31 Doppler velocities also agreed 
very well with the DWR-G displacements. The significant wave height had slightly 
higher error than that calculated from the GT-31 DGPS displacements, but was still 
within approximately 4% (0.14m) of the DWR-G buoy. The mean wave period differed 
by less than 4% (0.28sec) and the mean wave direction differed by less than 0.1°. Similar 
to the GT-31 DGPS displacements, the directional spread was within approximately 7% 
(~1.5°) of the DWR-G buoy. 


Table 11. Bulk wave statistics calculated from the two horizontal positioning methods 
are compared with the DWR-G results. 


Horizontal Positioning Comparison 


a x P x 
Positioning source : hide vd shalt Diff erat sented Diff 
Displacements | Displacements Velocities 





Mean wave period (sec) 8.4 8.9 
Directional spread (deg) a eee ee ee 





A spectral comparison of the two horizontal positioning techniques of the GT-31 
was also performed following the methods in Pearman et al. (2013). The wave energy 
spectra (top panel of Figure 28) were calculated using the horizontal surface motions 
from the DWR-G, the GT-31 DGPS displacements and the GT-31 Doppler velocities. All 
three estimates have similar spectral characteristics and spectral energy levels, 
particularly between 0.085Hz — 0.18Hz, which encompasses the typical swell-sea band. 
At frequencies below 0.085Hz (long period waves), both GT-31 estimates have higher 
spectral energy levels compared to the DWR-G. At frequencies above 0.18Hz (short 
period waves), the GT-31 DGPS displacements observe reduced spectral energy levels, 


while the GT-31 Doppler velocities observe similar spectral energy levels as the DWR-G. 


Spectral estimates of the mean wave direction calculated from 2nd-order Fourier 
moments are compared in the middle panel of Figure 28. All three estimates agree 
extremely well for frequencies between 0.078Hz — 0.24Hz. At frequencies less than 
0.078Hz, the DWR-G directional spectra has some erratic direction observations, which 


are likely due to the 180° ambiguity introduced by the 2nd-order Fourier moments as 
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noted by Herbers et al. (2012). At frequencies above 0.24Hz the GT-31 DGPS 
displacements show more variability than the DWR-G and the GT-31 Doppler velocities. 


Spectral estimates of the directional spread calculated from 2nd-order Fourier 
moments are compared in the bottom panel of Figure 28. The three estimates show 
similar variations with frequency, with the narrowest spread near the 0.1Hz swell peak 


frequency. At higher and lower frequencies, both GT-31 estimates are biased high by a 
few degrees relative to the DWR-G. 


Datawell Intra-Comparison 
Wave Energy and Directional Spectra from Horizontal Motions 
T T x T 








a 
Datawell Differential GPS 
—— GT-31 Differential GPS 
—— GT-31 Doppler Velocity 




















a 
a 


Horizontal Wave Spectrum (m@/Hz) 

















a 
D 
= 
c 
= 
Gb 
2 
5 
@ 
> 
2 
(= 
ig 
Cy 
zi 





f=) 
a 





Directional Spread (deg) 











Frequency (Hz) 





Figure 28. Wave spectra calculated from horizontal surface motions. Top panel: wave 
energy spectra. Middle panel: mean wave direction. Bottom panel: directional 
spreading. 


This confirms the results from Herbers et al. (2012) that the GT-31 DGPS 
displacements can provide similar bulk wave statistics and spectra as that of the DWR-G, 
particularly in the swell-sea frequency band. This analysis also shows that the statistics 
and spectra calculated from the GT-31 Doppler velocities are in excellent agreement with 


the DWR-G estimates. Further, the GT-31 Doppler velocities perform better than the GT- 
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31 DGPS displacements at higher frequencies, making this horizontal positioning 
technique more advantageous for observing the wave-induced surface motion across a 


wide frequency range. 


3. Vertical Wave Motion 


The sensor configuration developed for the WRD provided two different 
techniques for measuring the vertical surface motions of the waves: traditional DGPS 
displacements and vertical accelerations. Similar to the horizontal positioning, the 


performance of these two techniques was compared to that of the well-tested DWR-G. 


A plot of the vertical displacement time series from the three sensors is shown in 
Figure 29. All three estimates show the same phase evolution, but there are large 
differences in the amplitude of the waves. The DWR-G displacements and the X6—2 
acceleration-based estimates agree well with some tendency for the X6—2 accelerations to 
yield larger wave heights. The faster sampling rate of the accelerometer compared to the 
DWR-G buoy may contribute to these differences by better resoling the wave crests. On 
the other hand, the accelerometer-based wave heights may be biased high by uncorrected 
slope errors (see Pearman et al. 2013 for more discussion). Conversely, the GT-31 DPGS 
measures significantly lower wave amplitudes as compared to the DWR-G buoy 
(Pearman et al. 2013). This is likely due to the poor vertical resolution of the single 


frequency GPS receiver in the GT-31. 
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Figure 29. Vertical displacement time series. 


Statistics of the measured vertical surface motions are listed in Table 12 for 
comparison. Bulk wave statistics from the GT-31 DGPS displacements showed 
extremely poor agreement with the DWR-G displacements. The significant wave height 
is biased low almost 45% (1.5m) as compared to that calculated by the DWR-G buoy. 
Significant errors are also noted in the mean wave period (10%, 0.83sec) and directional 
spread (51%, 13.7°). However, the mean wave direction differs by only 1.3° from the 
direction measured by the DWR-G buoy, indicating that the phase evolution of the wave 


field is still well resolved. 


The bulk wave statistics from the X6—2 accelerations are in much better 
agreement with the DWR-G buoy. The significant wave heights differed by less than 6% 
(0.19m) and the mean wave period and the mean wave direction also agree well. On the 
other hand, the directional spread is almost 28% (7.42°) higher than the DWR-G buoy 


estimate, suggesting higher noise levels in the accelerometer measurements. 
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Table 12. Bulk wave statistics calculated from the three vertical positioning methods 
with differences compared to the DWR-G displacements. 


Vertical Positioning Comparison 


Positionin : DWR-G GT-31 DGPS 
ils cist Displacements | Displacements 


Significant wave height (m) 44.37 % 


Mean wave period (sec) | ss | tt | 9% | 
Mean wave direction (deg) | 309] 0 | 30" 
Directional spread (deg) 





A spectral comparison of the two vertical positioning techniques was also 
performed following the methods in Pearman et al. (2013). The wave energy spectra (top 
panel of Figure 30) were calculated using the vertical surface motions from the DWR-G, 
GT-31 DGPS displacements and X6—2 accelerations. All three positioning techniques 
have similar spectral characteristics; however they differ in their spectral energy levels. 
The DWR-G and the X6—2 accelerations agree very well for frequencies between 0.08Hz 
— 0.18Hz, but diverge slightly at frequencies greater than 0.18Hz (short period waves), 
where the X6-—2 accelerations yield a slightly more energetic spectrum, possibly because 
the X6—2 accelerometer is more sensitive to nonlinear (i.e. near breaking) high frequency 
waves. Conversely, the GT-31 DGPS displacements show much lower spectral energy 
levels at all frequencies in the wind wave and swell band (Herbers et al. 2012 and 


Pearman et al. 2013). 


Spectral estimates of the mean wave direction calculated from Ist-order Fourier 
moments are compared in the middle panel of Figure 30. Generally, all three estimates of 
the mean wave direction agree well. At frequencies less than 0.075Hz, where energy 
levels are relatively low, there are large discrepancies. Neither of the WRD vertical 
measurement techniques accurately compare to the DWR-G at these low frequencies. At 
frequencies greater than 0.22Hz, the X6—2 accelerations agree well with the DWR-G 


displacements, while the GT-31 DGPS displacements show some variability in direction. 


Spectral estimates of the directional spread calculated from Ist-order Fourier 
moments are compared in the bottom panel of Figure 30. The directional spread shows 
large variability in the degree of spread among the three positioning techniques. The GT- 
31 DGPS displacements and X6—2 accelerations observe consistently larger spreading 
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across the entire frequency band as compared to the DWR-G displacements. However, 
the X6—2 acceleration spectrum does agree better than the GT-31 DGPS displacement 


spectrum, particularly for frequencies greater than 0.08Hz. 
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Figure 30. Wave spectra calculated from vertical surface motions. Top panel: wave 
energy spectra. Middle panel: mean wave direction spectra. Bottom panel: 
directional spreading. 


This analysis shows that the bulk wave statistics and wave energy spectra 
calculated using the GT-31 DGPS vertical displacements do not compare well to those 
calculated by the DWR-G buoy, which is consistent with previous results from Herbers et 
al. (2012). However, these results show that the X6—2 accelerometer does provide 
reliable observations of the vertical wave-induced surface motion, which can be used in 
conjunction with the GT-31 horizontal measurements to obtain bulk wave statistics and 


wave spectra. 


Overall, the combined use of the horizontal Doppler velocities from the GT-31 


GPS receiver and the vertical accelerations from the X6—2 accelerometer provide an 


ny) 


accurate and reliable sensor configuration for observing the wave-induced fluid motions. 
Collocated on the same buoy, these sensors accurately observe the surface motions in the 
swell and wind-wave frequency band, producing robust wave and spectral statistics that 


are similar to those provided by widely used Datawell buoys. 


C. INTER-BUOY COMPARISON: DWR-G VS. WRD 


The second step of the validation involved a statistical and spectral comparison of 
the WRD buoy design. Three WRD buoys and one DWR-G buoy were outfitted with the 
same sensor configuration and were deployed in the same homogenous deep-water wave 


field. 


The sensor validation showed that the most reliable statistics for the swell-sea and 
wind-wave frequency band were obtained using the GT-31 horizontal Doppler velocities 
in conjunction with X6—2 vertical accelerations. The three WRD buoys (A06, AO8 and 
A17) were processed using these data streams. In addition, the GT-31 and X6—2 sensor 
configuration onboard the DWR-G was also included in the analysis and processed the 


same way as the stand-alone WRD buoys. 


Statistics calculated from the DWR-G and WRD buoys using the same sensors 
(GT-31 and X6-2) are listed in Table 13 for comparison. The bulk wave statistics from all 
three WRD buoys agreed very well with the statistics calculated from the DWR-G. The 
significant wave height calculated by the WRD buoys differed by an average of less than 
7% (0.24m) when compared to the DWR-G buoy. Other wave parameters show similar 
good agreement with differences of 2% (0.13sec) for the mean wave period, less than 1° 
for the mean wave direction, and 6% (1.5°) for the directional spread, as compared to the 


DWR-G buoy. 


The bulk wave statistics showed little difference between those calculated from 
vertical surface motions and those calculated from horizontal surface motions. The 
average significant wave height of the WRD buoys calculated from vertical motions was 
within 6% (0.21m) of the DWR-G buoy while that calculated from horizontal motions 


was within 7.5% (0.27m). Similarly, the average directional spread of the WRD buoys 
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calculated from vertical motions was within 5% (1.6°) of the DWR-G buoy while the 


directional spread calculated from horizontal motions was within 6% (1.4°). 


Table 13. Bulk wave statistics calculated from the three WRD buoys and DWR-G 
buoy using GT-31 Doppler velocities and X6—2 accelerations with 
differences compared to the DWR-G buoy. 


WRD COMPARISON to DWR-G (Doppler velocity/Accelerations) 
[_veriatonietictons Tpenwet [ wenves [om | weosae [et _[ weeasy [oat [ wenaeerees [om 
Significant wave height (m) i 
Mean wave period (sec) 
Mean wave direction (deg) 
Directional spread (deg) 











WRD cariga pon 
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Mean wave direction (deg) 

















Directional spread (deg) 


A further comparison of the bulk wave statistics showed that the average statistics 
calculated from the WRD buoys compared very well to the statistics calculated from the 
DWR-G displacements (see Table 14). The average significant wave height of the WRD 
buoys was with 12% (0.40m) when compared to the significant wave height from the 
DWR-G displacements. The larger positive bias (compared with Table 13) shows that 
both the buoy response and the sensor characteristics contribute to these differences. The 
average wave period, wave direction and directional spread of the WRD buoys are within 
1% (0.16sec), 1° and 6° respectively, when compared to the statistics from the DWR-G 


displacements. 


Table 14. | Bulk wave statistics calculated from the three WRD buoys using GT-31 
Doppler velocities and X6—2 accelerations with differences compared to the 
DWR-G DGPS displacements. 


WRD COMPARISON to DWR-G (Displacements) 


[__vertenloroal eters [ Deewot [ wan-aee [om _[ wees | per won-at7 | oar_T wan Areree [oat] 
Mean wave period (sec) 8: 2 2.09 % 8.2 -0.41% 8.2 -0.06 % 8.3 0.54 % 
Directional spread (deg) a 5.29° 34 7.02° 32 5.05° 33 5.79° 








Horizontal Orbital Motions | Datawell | WRD-AQ6 Diff wrb-Aos | Diff | WRD-A17 viff | WRDAverage | Diff | 


Significant wave height (m) 3.43 3.34 -2.69 % 3.19 -7.16 % 3.39 -1.28 % 3.31 -3.71% 
Mean wave period (sec) 8.4 8.7 3.59% 8.6 2.40 % 8.7 2.96 % 8.7 2.99 % 
Mean wave direction (deg) 310 312 1.79° 311 0.72° 310 -0.18° 311 0.77° 


Directional spread (dee) 21 013° [2 | 079" | 20 0.36" 
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A spectral comparison of the three WRD buoys and the DWR-G buoy was also 
performed following the methods in Pearman et al. (2013). Spectra calculated from 
vertical surface motions are shown in Figure 30. The vertical wave energy spectrum (top 
panel of Figure 31) was calculated using the surface motions from the DWR-G 
displacements and the X6-2 accelerations on the WRD buoys. The DWR-G, WRD-A06, 
WRD-AO8 and WRD-A17 buoys show excellent agreement in spectral peaks, 
characteristics and spectral energy levels for frequencies greater than 0.08Hz. At 
frequencies less than 0.08Hz, the WRD buoy spectra show similar spectral energy levels 


that are higher than the DWR-G estimates. 


Spectral estimates of the mean wave direction calculated from 1st-order Fourier 
moments are compared in the middle panel of Figure 31. Again, the WRD buoys and 
DWR-G buoy show excellent agreement for all frequencies greater than 0.075Hz. At 
frequencies less than 0.075Hz, the directional spectra for all four buoys diverge 
indicating that this sensor and buoy configuration cannot accurately resolve these long- 


period swell directions. 


Finally, spectral estimates of the directional spreading calculated from Ist-order 
Fourier moments are compared in the bottom panel of Figure 31. Again, the WRD buoys 
and DWR-G buoy spectra show excellent agreement for all frequencies greater than 
0.065Hz. At frequencies less than 0.065Hz, the WRD buoys agree well with each other, 
but display slightly higher degrees of spreading than the DWR-G. 
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Figure 31. Wave spectra calculated from vertical surface motions. Top panel: wave 
energy spectra. Middle panel: mean wave direction spectra. Bottom panel: 
directional spreading. 


Spectra calculated from horizontal surface motions are shown in Figure 32. The 
horizontal wave energy spectra (top panel of Figure 32) were calculated using the 
horizontal DWR-G displacements and GT-31 Doppler velocities, applying the linear 
transfer functions from Table 6. Similar to the vertical spectra, the DWR-G buoy, WRD- 
A06, WRD-A08 and WRD-A17 show excellent agreement in spectral peaks, 
characteristics and spectral energy levels for the entire frequency band. There is some 
variability at frequencies less than 0.07Hz and greater than 0.25Hz where the DWR-G 


spectrum has slightly higher spectral energy levels. 


Spectral estimates of the mean wave direction calculated from 2nd-order Fourier 
moments are compared in the middle panel of Figure 32. Again, the WRD buoys and 
DWR-G buoy show excellent agreement for frequencies between 0.08Hz — 0.2Hz. Unlike 
the vertical motions, the horizontal motions show greatly improved agreement at low 
frequencies (less than 0.08Hz); however, they also show larger variability at frequencies 


greater than 0.2Hz. 
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Finally, spectral estimates of the directional spreading calculated from 2nd-order 
Fourier moments are compared in the bottom panel of Figure 32. Again, the WRD buoys 
and DWR-G buoy spectra show good agreement in their spectral characteristics across 


the entire frequency band; however, there is increased variability at frequencies greater 
0.085Hz. 
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Figure 32. Wave spectra calculated from horizontal surface displacements. Top panel: 
wave frequency spectra. Middle panel: mean wave direction spectra. Bottom 
panel: directional spreading. 


The statistical and spectral comparison confirms that the WRD buoy and sensor 
configuration provides similar bulk wave statistics and wave spectra as the widely used 
DWR-G buoys, particularly in the dominant swell and wind-wave frequency band (0.07- 
0.25Hz). Table 15 summarizes the WRD system performance compared to the DWR-G. 
The total system performance is comprised of errors/bias introduced by the sensor 
configuration (evaluated in Section 4.A) and the buoy design (evaluated in Section 4.B). 
The overall system performance is slightly better using horizontal surface motions; 


however, the vertical surface motions still provide adequate statistics as well. The 
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significant wave height calculated from the horizontal motions is within 3.7% of the 
DWR-G while the wave height calculated from the vertical motions is within 12% of the 
DWR-G. 


Table 15. | WRD system performance. 


DWR-G (GPS) vs. DWR-G (Doppler Velocity / DWR-G (GPS) vs. 
Sensor Configuration DWR-G (Doppler Velocity Acceleration) vs. WRD (Doppler 
/ Acceleration) WRD (Doppler Velocity / Velocity / 
Acceleration) Acceleration) 


Vertical Motions 


Significant wave height (m) 11.94 % 


Mean wave period (sec) ec ee 2.75 % 0.54 % 
Meanwavedirection (deg) | one | aasé 0.97" 


Horizontal Motions 

|Meanwaveperiod(sec) | = 3.32% =| 032% 2.99 % 
0.77" 
0.10" 


Directional spread (deg) 42° a 5.79° 


-7.51% -3.71% 
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V. APPLICATION OF WRD BUOY IN THE NEARSHORE 


The WRD buoys were validated to provide similar bulk wave statistics and 
spectra as the DWR-G in a homogenous offshore wave field. In order to examine their 
applicability in the nearshore environment, an array of WRD buoys was deployed on four 
cruises within the Golden Gate channel and the mouth of San Francisco Bay. These 
buoys drifted with the surface tidal current and observed the spatial and temporal 
evolution of the wave field across this dynamic nearshore environment. Two of the four 


nearshore deployments are discussed here. 


A. DRIFTER TRACKS 


The first nearshore deployment was conducted on February 15, 2012. 20 WRD 
buoys were deployed from underneath the Golden Gate Bridge on an ebb-tidal current 
and allowed to drift out of the channel with the current; however only 13 WRD buoys 
were recovered at various times over the next 24 hours. WRD deployment and retrieval 
times are listed in Table 16 along with the applicable tidal currents that each buoy 


experienced. 
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Table 16. February 15 WRD buoy deployment times. 


February 15, 2012 Drifter Deployment 


= —— Tidal Current 


Time Time 


Ebb #1 | Flood #1 


Ebb #1 | Flood #1 


wRD-Ai3 | 2/15 | 17: ae 
WRD-A14 07:53 | Ebb #i | Flood # | 
WRD-A15 | 2/15 2/15 | 23:21 


WRD-A17 | 2/15 2/15 rebb#i| x | 
WRD-A18 : aa oe 


WRD-A19 | 2/15 | 18:10 | 2/16 Ebb #1 | Flood #1 


WRD-A21 | 2/15 





The drift tracks for this deployment are shown in Figure 33. Since the 
predominant wind was out of the NW (318°) at 10m/s, the drift tracks demonstrate that 
the buoys were well coupled to the surface current and not strongly affected by the wind. 
This deployment demonstrates the possibility of utilizing a large array of these 
Lagrangian drifters in order to observe the spatial evolution of the wave field in the 
nearshore environment, particularly the interaction with bathymetry and variable surface 


tidal currents. 
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Figure 33. Buoy drift tracks from the February 15 deployment indicating the tidal stage 
on the drift tracks. 


WRD-AO1, WRD-A02, WRD-A03 and WRD-A06 were deployed first and 
experienced the strongest ebb current. These buoys drifted out of the channel and well 
past the ebb-tidal shoal before the ebb current abated. Offshore of the ebb-tidal shoal, the 
currents were weak, and all four of these buoys remained offshore, oscillating in 
successive flood and ebb currents. However, none of these currents were strong enough 


to bring the buoys back toward shore. 


Conversely, WRD-A14, WRD-A16 and WRD-A19 were deployed later in the 
day, on a weaker ebb current. These buoys only drifted out as far as the ebb-tidal shoal 
before the tide shifted. The subsequent flood current near the shoal was strong enough to 
affect the buoys, bringing them back toward the channel before WRD-A16 and WRD- 


A19 were retrieved. 


WRD-A14 continued to drift toward shore on the flood current. This buoy was 
close enough to shore during the next tide change, that during this ebb current, the buoy 


began drifting in an alongshore current south of the channel. This southward flowing 
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alongshore current is also seen in the drift track of WRD-A17 which drifted south as it 
exited the channel on the initial ebb current. These two buoy drift tracks identify an ebb- 


tidally driven alongshore current south of the channel discussed in Shi et al. (2011). 


On the north side of the channel, WRD-AO7 drifted north on the initial ebb 
current as it exited the channel and drifted out over the ebb-tidal shoal. As the tide 
changed, the buoy drifted back toward shore on the flood current. It continued to drift 
southeastward with the flood current moving in and out of the surf zone, until it drifted 
offshore with the next ebb current. The buoy drifted offshore and northwestward during 
the ebb current. This buoy drift track identifies a tidally-driven eddy circulation that 


forms over the northern reach of the ebb-tidal shoal outside of the channel. 


B. SPATIAL VARIABILITY ANALYSIS 
1. February 15, 2012 


The data from this deployment was analyzed for the period of the first ebb current 
(15:14 UTC — 21:36 UTC). The buoy drift tracks show the spatial coverage of the WRD 
array during this time (see Figure 34). In order to study some of the dynamic features of 
the nearshore environment, five distinct spatial regions of the dataset were identified: 
offshore of the ebb-tidal shoal, over the ebb-tidal shoal, inshore of the ebb-tidal shoal, in 
the vicinity of Mile Rock (the entrance to the Golden Gate channel) and inside the 
channel (see Figure 34). These regions were identified for having similar wave statistics. 
The bulk wave statistics for these regions were averaged for all of the available WRD 


buoys and are listed in Table 17. 
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Figure 34. Buoy drift tracks from the February 15 deployment. The left panel shows the 
overall spatial distribution. The right panel indicates the five spatial regions. 


Table 17. Bulk wave statistics averaged over all of the WRD buoys for the five spatial 
regions of the February 15 deployment. 
WRD Golden Gate Drift Deployment — Feb 15, 2012 


Offshore Shoal Inshore Channel Entrance Channel 
Std Dev Mean Std Dev Mean Std Dev Mean Std Dev Mean Std Dev 


Vertical Orbital Motions 





Depth (m) 
Mean drift rate (m/s) 





Significant wave height (m) : : b 0.24 


Mean wave period (sec) d : 5 











The buoys drifted out the channel, in the opposite direction of the approaching 
wave field. The offshore wave field (measured by NDBC station 46026) is characterized 
by a significant wave height of approximately 2.7m coming out of the WNW (286°) 
during this deployment period (from Table 4). 


As the wave field propagates from offshore to over the ebb-tidal shoal, the 
average water depth decreases significantly from 27m to 12m. The tidal current offshore 
of the shoal is relatively weak as evidenced by a low buoy drift rate of 0.5m/s (1.0kts). 
The significant wave height in this region increases by 0.27m. This is qualitatively 
consistent with the linear shoaling effect resulting in an increase of wave energy due to 
the decreasing group speed over the shoal. Additionally, the mean wave direction in this 


region changes slightly by coming from 280° to coming from 274°. This appears 
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consistent with refraction over the bathymetry as the wave field propagates over the 


dredged channel and onto the ebb-tidal shoal. 


As the wave field propagates past the shoal, the average water depth increases 
from 12m to 20m and the ebb-tidal current steadily increases, reaching an average speed 
of 2.3kts (using the buoy drift rate). The significant wave height decreases by 0.57m in 
this region which is consistent with bathymetric effects of the shoal. Much of the swell 
energy is likely dissipated over the shoal (Pearman et al. 2013) and the shoaling effect 
reverses as the wave field propagates towards deeper water. The mean wave direction 
only slightly changes direction back to coming from 276°. The directional spreading in 
this region also increases by 8° in this region. Since the bathymetry of the shoal continues 
to the southeast, these two directional effects suggest that the wave field over the 
southern reach of the shoal is beginning to refract back toward the channel and onto the 


opposing ebb-tidal current. 


As the wave field propagates toward the channel entrance, the average water 
depth continues to increase (to 40m) while the ebb-tidal current also continues to increase 
in strength, reaching an average speed of 2.9kts. The significant wave height in this 
region decreases by 0.38m while the mean wave direction and directional spreading 
remain constant from the inshore region. There is no obvious signature of any expected 
wave-current interactions, and the temporal and spatial variability of the tidal currents 
complicate the analysis in this region. The buoy observations occurred during peak ebb 
currents (early in the deployment), and very weak ebb currents (toward the end of the 
deployment). Additionally, the buoy drift tracks encompassed nearly the entire width of 
the channel, capturing observations inside the tidal jet, and the northern and southern 


boundaries of the jet. 


Finally, as the wave field propagates into the center of the channel, the average 
water depth increases significantly (to 92m) and the ebb current reaches its strongest 
speed, increasing to 3.1kts (from buoy drift rate). The significant wave height again 
decreases by 0.36m, which is consistent with the decrease in wave energy as the wave 
field spreads out, traveling faster in the deeper water of the channel. Inside the channel, 


the mean wave direction changes nearly 20° to the southeast (coming from 295°) as the 
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waves begin to exhibit some refractive effects from the nearby coastline. The directional 
spreading of the wave field also increases by 6° as the waves react to the refractive 
effects of the coastline as well as some refractive effects from the ebb current, which is 


strongest for all of the buoys in this region. 


2. April 27, 2012 


The second nearshore deployment was conducted on April 27, 2012. Similar to 
the February deployment, a large array of 29 WRD buoys was deployed from underneath 
the Golden Gate Bridge. However, during this deployment, the buoys were deployed in 
the center of the tidal jet as measured by the ADCP onboard R/V Questuary. The buoy 
drift tracks show the spatial coverage of the WRD array (see Figure 34). This deployment 
provided a unique synoptic view of the wave field as it progressed from offshore. The 
data was analyzed for the same five spatial regions used for the February dataset (see 
Figure 35). The bulk wave statistics for these regions were averaged for all of the 


applicable WRD buoys and are listed in Table 18. 
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Figure 35. Buoy drift tracks from the April 27 deployment. The left panel shows the 
overall spatial distribution. The right panel indicates the five spatial regions. 
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Table 18. Bulk wave statistics averaged over all of the WRD buoys for the five spatial 
regions of the April 27 deployment. 


WRD Golden Gate Drift Deployment — Apr 27, 2012 
Offshore Shoal Inshore Channel Entrance Channel 
Std Dev Mean Std Dev Mean Std Dev Mean Std Dev Mean Std Dev 


Depth (m) 28 12 19 36 
Mean drift rate (m/s) 0.7 1.4 1.8 1.9 


Vertical Orbital Motions 











Significant wave height (m) A L 2.80 7 2.14 0.51 
Mean wave period (sec) 
Mean wave direction (deg) 

















The offshore wave field measured by NDBC Station 46042 was similar to 
February 15, 2012 with a significant wave height of 2.6m and a wave field coming out of 
the WNW (288°) (from Table 5). Again, the buoys were well coupled to the surface 
current as the buoys drifted straight out of the channel despite having prevailing surface 


winds of 11m/s out of the northwest (315°). 


As the wave field propagates from offshore to over the shoal, the average water 
depth decreases significantly from 28m to 12m. Again, the offshore current is relatively 
weak as evidenced by a buoy drift rate of 0.7m/s (1.4kts). Similar to the February 
deployment, the significant wave height increases by 0.29m as the wave field propagates 
over the shoal as a result of the bathymetric interaction. The mean wave direction 
changes slightly by turning nearly 4° to the east (coming from 273°) and the directional 
spreading narrows by more than 4°. These directional changes are qualitatively consistent 


with refraction over the shoal. 


As the wave field propagates off the shoal, the average water depth increases from 
12m to 19m. The effects of the ebb-tidal current begin to increase as the buoy drift rate 
increases to 1.8m/s (3.5kts). Similar to February, the significant wave height decreases by 
0.66m as much of the swell energy is dissipated from the bathymetric interaction of the 
shoal (Pearman et al. 2013). The decrease in swell energy is also observed in a decrease 
in the mean wave period from 8.5sec to 7sec indicating that the wave field has 
diminished contributions from long-period swell. As seen in February, the wave 
directions do not continue to refract to the southeast along the shoal; instead, the mean 
wave direction remains eastward (coming from 273°). This lack of directional change 
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accompanied with a 9° increase in the directional spread indicates influence of the 
Pp Pp: 


opposing ebb-tidal current resulting in the wave field refracting back toward the current. 


As the wave field propagates toward the channel entrance, the water depth 
continues to increase (to 36m), while the ebb current reaches its maximum speed of 
1.9m/s (3.7kts) indicating that the ebb-tidal current is strongest in this region. The mean 
wave direction changes by 10° (coming from 283°) as the wave field observed by the 
buoys refracts into the channel. This results in the wave field propagating against the 
strong, ebb current. The effects of the ebb current are also exhibited in the 4° narrowing 


of the directional spreading, indicating that the wave field is refracting in on the current. 


Finally, as the wave field propagates into the center of the channel, the water 
depth increases significantly (to 97m) and the buoy drift rate abates slightly to 1.6m/s 
(3.1kts). Much of the wave energy has either been dissipated or blocked by the opposing 
ebb current at the channel entrance, which results in a 0.63m decrease in the significant 
wave. The mean wave direction continues to change by nearly 20° (coming from 305°) 


inside the channel, possibly owing to sheltering effects of the nearby coastline. 
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VI. CONCLUSIONS 


The nearshore environment, particularly around a narrow tidal inlet, provides 
many challenges for observing ocean waves with traditional sensors. Migrating sand 
waves and sediment transport provide a continual shifting bathymetry that is not 
conducive for moorings of traditional Eulerian wave buoys, while strong (tidal) currents 
can strain surface buoy mooring lines. Additionally, many of these tidal inlets and 
harbors, such as San Francisco Bay, are busy shipping ports with a steady flow of large 


commercial shipping. 


These same challenges that make wave measurements difficult in inlets also result 
in strongly variable wave statistics, understanding of which is important for circulation 
and transport processes in the inlet, and navigation. A Lagrangian GPS drifter buoy is 
well-suited for the measurement of the wave field in these environments, but a modern 
Datawell Waverider-GPS buoy can cost as much as $14,000, making it prohibitively 
expensive to deploy in large numbers in an environment where one might be lost or 


damaged. 


In this study, an inexpensive (less than $1,000) Wave Resolving Drifter (WRD) 
buoy developed at NPS (Herbers et al., 2012) has been improved through the use of 
horizontal GPS Doppler velocity measurements and the addition of a vertical-axis 
accelerometer. Transfer functions relating these various measurements to sea surface 
elevations (wave height) were applied and used to estimate the wave frequency and 
directional spectra observed by the buoy. These spectral estimates were then used to 


calculate the bulk wave statistics of the wave field. 


The Doppler velocities from the GT-31 GPS receiver provide an order of 
magnitude improvement in the horizontal measurement over the traditional DGPS 
displacement measurements. An analysis was performed comparing the wave spectra and 
statistics from the traditional DGPS displacement measurements from the GT-31 and 
Doppler velocity measurements from the GT-31 to those of the well-tested Datawell 


DWR-G buoy. Although both GT-31 measurements provided similar statistics as the 
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DWER-G in the low frequency swell band; the Doppler velocity measurements agreed 
much better than the DGPS displacements with the Datawell estimates in the higher 


frequency wind wave band. 


The addition of the vertical accelerometer to the WRD buoy design provided a 
significant improvement in the accuracy of vertical observations, which was previously 
lacking due to the poor vertical resolution of traditional DGPS displacement 
measurements. An analysis of the wave spectra and statistics demonstrates excellent 
agreement between the accelerometer and the DWR-G across the entire frequency band, 
while the traditional DGPS displacement measurements (from the GT-31) performed 


poorly as compared to the DWR-G. 


This improved accuracy of the vertical observations allows for greater flexibility 
in calculating the bulk wave statistics. Previous work relied on horizontal measurements 
only, which yielded accurate results, but contain a 180° ambiguity in the mean wave 
direction (Herbers et al. 2012). The improved vertical accuracy from the accelerometer 
allows for the reliable use of vertical measurements as well, which removes this 180° 


directional ambiguity. 


The validation of the WRD buoy, performed in deep water, free of any 
bathymetric and coastal sheltering effects, demonstrates its capability to collect robust 
estimates of wave spectra and statistics comparable to the industry-standard Datawell 


Waverider GPS buoy (DWR-G). 


A large array of WRD buoys (20-30 buoys) was deployed during four cruises in 
and around the San Francisco Bay entrance. The buoys were well coupled to the surface 
current and allowed to drift out of the Golden Gate channel on the outgoing ebb current. 
The buoys drifted through a region of strong wave-current interactions and then over an 
ebb-tidal shoal, which resulted strong bathymetric effects (shoaling and refraction). 
Additionally, the drift tracks identify tidally-induced eddy circulations and alongshore 
currents near the mouth of the channel (the exit region of the tidal inlet). The spatial 


evolution of the wave field was analyzed on the ebb current of two deployments by 
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separating the time series geographically in order to statistically examine the wave field 


in dynamically different regions of the nearshore environment. 


As the wave field propagated from offshore toward the ebb-tidal shoal, the WRD 
buoys observed an increase in significant wave height associated with the bathymetric 
interaction with the ebb-tidal shoal. The buoys also measured the distinct decrease in 
spectral energy in the swell band over the shoal that is likely caused by bottom friction 
(Pearman et al. 2013). As the wave field propagated off the shoal and into the inshore 
region, the buoys observed a decrease in significant wave height due to the dispersion of 
wave energy in the deeper water. However, they also observed some refractive effects 
and increased directional spreading of the wave field associated with the influence of the 
opposing ebb current. As the waves propagated into the channel, the buoys observed a 
significant decrease in significant wave height that is likely caused by a combination of 
the current blocking high frequency wind waves and the adjacent coastline sheltering the 


swell arrivals. 


This study has shown that an inexpensive drifting buoy equipped with GPS and 
an accelerometer can provide reliable wave statistics. An array of these buoys can be 
used to observe the spatial and temporal evolution of a wave field in a dynamic nearshore 
environment. Analysis of these datasets is currently underway to study wave-current and 
wave-bottom interactions. The wave statistics, drift tracks and current velocities are also 
intended to be compared to the output of a coastal wave and circulation model for 
comparison. Future research is planned for deploying large arrays of these WRD buoys 
near the mouth of the Columbia River in Oregon later this year in order to further study 


wave-current interactions in the presence of tidally-enhanced river outflow. 
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